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Abstract. We discuss the thermoelectric (TE) transport in billiard systems of interacting particles, 
coupled to stochastic particle reservoirs. Recently in [1], analytical exact expressions for the TE 
transport of noninteracting gases of polyatomic molecules were obtained, and a novel microscopic 
mechanism for the increase of thermoelectric efficiency described. After briefly reviewing the 
derivation of [1], in this paper we focus on the effects that the particle-particle interaction has on the 
TE efficiency. We show that interaction reduces the maximal thermodynamic efficiency. However, 
the mechanism for the efficiency's increase towards its Carnot upper limit, described in [1], remains 
unaffected. 
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1. INTRODUCTION 

Thennoelectricity concerns the conversion of temperature differences into electric po- 
tential or vice-versa. It can be used to perform useful electrical work or to pump heat 
from cold to hot place, thus performing refrigeration. Although thermoelectricity was 
discovered almost 200 years ago, a strong interest of the scientific community arose only 
in the 1950's when Abram loffe discovered that doped semiconductors exhibit relatively 
large thermoelectric effect. This initiated an intense research activity in semiconductors 
physics which was not motivated by microelectronics but by the loffe suggestion that 
home refrigerators could be built with semiconductors [2, 3]. As a result of these efforts 
the thermoelectric material Bi2Te3 was developed for commercial purposes. However 
this activity lasted only few years until the mid 1960's since, in spite of all efforts and 
consideration of all type of semiconductors, it turned out that thermoelectric refrigera- 
tors have still poor efficiency as compared to compressor based refrigerators. Nowadays 
Peltier refrigerators are mainly used in situations in which reliability and quiet opera- 
tion, and not the cost and conversion efficiency, is the main concern, like equipments 
in medical applications, space probes etc. In the last decade there has been an increas- 
ing pressure to find better thermoelectric materials with higher efficiency. The reason is 
the strong environmental concern about chlorofluorocarbons used in most compressor- 
based refrigerators. Also the possibility to generate electric power from waste heat using 
thermoelectric effect is becoming more and more interesting [4, 2, 3]. 

The suitability of a thermoelectric material for energy conversion or electronic refrig- 
eration is evaluated by the thermoelectric figure of merit Z, 



where a is the coefficient of electric conductivity, S is the Seebeck coefficient and 
K is the thermal conductivity. The Seebeck coefficient S, also called thermopower, 
is a measure of the magnitude of an induced thermoelectric voltage in response to a 
temperature difference across the material. 

For a given material, and a pair of temperatures Tn and Tq of hot and cold thermal 
baths respectively, Z is related to the efficiency r] of converting the heat current Jq (be- 
tween the baths) into the electric power P which is generated by attaching a thermoelec- 
tric element to an optimal Ohmic impedance. Namely, in the linear regime: 
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where r/camot = 1 — 7c/7h is the Camot efficiency and T = (Th + Tc)/2. Thus a good 
thermoelectric device is characterized by a large value of the non-dimensional figure of 
merit ZT. 

Since the 1960's many materials have been investigated but the maximum value 
found for ZT was achieved for the (Bii_xSbx)2(Sei-jTey)3 alloy family with ZT ^ 1. 
However, values ZT > 3 are considered to be essential for thermoelectrics to compete 
in efficiency with mechanical power generation and refrigeration at room temperatures. 
The efforts recently focused on a bulk of new advanced thermoelectric materials and on 
low-dimensional materials, and only a small increment of the efficiency, ZT < 2.6, has 
been obtained [4]. 

One of the main reasons for this partial success is a limited understanding of the 
possible microscopic mechanisms leading to the increase of ZT, with few exceptions 
[5]. From a dynamical point of view, cross effects in transport have been barely studied 
[6, 7]. So far, the challenge lies in engineering a material for which the values of S, o 
and Jc can be controlled independently. However, the different transport coefficients are 
interdependent, making optimization extremely difficult. 

In a recent paper, we have studied the thermoelectric process in a gas of non- 
interacting polyatomic molecules [1]. We showed that large values of ZT, in principle 
approaching to Camot' s efficiency, are obtained when the number of the molecule's in- 
ternal degrees of freedom (d.o.f.) is increased. Using the rotating-disk interaction intro- 
duced in [7], in this paper, we study the effect that a generic particle-particle interaction 
has on the TE efficiency of these systems. In the following section 2, we review briefly 
the microscopic expressions for the TE process in ergodic ideal gases obtained in [1]. In 
the section 3, we study numerically, the modifications to these expressions, due to the 
interaction. Our conclusions appear in section 4. 

2. THERMOELECTRIC TRANSPORT 

In the linear response regime (see e.g. [8]), the transport equations for a thermoelectric 
material give the heat current Jq and the electric current Jg through an homogeneous 
sample subjected to a temperature gradient dxT and a electrochemical potential gradient 



dxlJ- as 



Jq= -id d^T-ToSd^fi , 
Jg = — oS dxT — o dxjl ■ 

Here and in what follows, we assume that the transport occurs along the x-direction and 
the temperature is given in units where the Boltzmann constant = I- 

The electrochemical potential is the sum of a chemical and an electric part /i = /i 
where jl is the chemical potential of the particles and, if e is the particle's charge, jig — 
is the work done by the particles against an external electric field S = — (9^0. From (3) 
the usual phenomenological relations follow: if the thermal gradient vanishes, (9^7 = 0, 
then Je = — cr d^^ = , since for an isothermal homogeneous system jU is uniform. 
If the electric current vanishes, — 0, then dx\i — S dxT, which is the definition of 
the Seebeck coefficient, and Jq = —KdxT where Jc = Jc' — TaS^ is the usual thermal 
conductivity (see e.g. [10]). 

From the theory of irreversible thermodynamics, /i and /ig cannot be determined 
separately; only their combination in p. appears in (3) [9]. Based on this equivalence, 
in what follows we take into account the chemical part only, i.e., p, = ji. 

Since in the linear regime in which dissipative effects, such as, e.g.. Joule heating, 
can be neglected and at low particle densities, thermodiffusion and thermoelectricity are 
equivalent process, we can study TE in terms of the thermodiffusion transport equations 



Ju — LuU dx ^ J' ^ J, 



1 ^^'^ 

where Ju and Jp = J^je are the energy and particle density currents, and Lj,p = Lp^ 
follows from the Onsager reciprocity relations. From the entropy balance equation for 
open systems [11] 

Ju = JQ + fiJp, (5) 

and substituting Jq in (3) in favor of Ju and comparing the resulting equations with (4) it 
follows that the TE transport coefficients can be written in terms of the thermodiffusion 
L-coefficients as 



e 



2- 1 detL „ 1 /L 



(^=^Lno, K=—^- , S= — {-—^-LL]. (6) 
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Furthermore, from Eq. (1), we obtain for the figure of merit 

~ detL • ^'^ 

Note that in Eqs. (6) and (7), T and jU are taken as the mean values. 

Consider now an ergodic gas of non-interacting, electrically neutral particles of mass 
m with Jint internal d.o.f. (rotational or vibrational), enclosed in a J dimensional con- 
tainer. To study the non-equilibrium state of such dilute poly-atomic gas we consider a 




FIGURE 1. The open Lorentz gas system and a typical particle's trajectory. The particles are scattered 
from fixed disks of radius R disposed in a triangular lattice at critical horizon, i.e., the width and height 
of the cells are Ax = 2R and Ay — 2W respectively, where W = 4-R/^/3 is the separation between the 
centers of the disks. The channel is coupled at the left and right boundaries to two thermochemical baths 
at temperatures Tl and Tr and chemical potentials /Xl and /Xr, respectively. 

chaotic billiard channel (like the one shown in Fig. 1) connected through openings of 
length A to two reservoirs of particles which are idealized as infinite chambers with the 
same poly-atomic gas at equilibrium density p and temperature T. From the reservoirs, 
particles are injected into the channel at a rate 7, as explained in the appendix A 

The particle injection rate 7 is related to the value of the chemical potential jl at the 
reservoirs which, for a gas of polyatomic molecules with a total of D = d + Jint d.o.f., at 
equilibrium density p and temperature T reads 

M = Mo + nn(^)=M^ + rin(^;(^) , (8) 

where /io and /ig are reference values of the chemical potential and the second equality 
is simply obtained after substitution of p from Eq. (17). Furthermore, energy is injected 
from the reservoirs at a rate £ = yT{D+1)/2 (see appendix A). 

Calling pt{l) the transmission probability of the channel of length /, the density cur- 
rents 7m, 7p for noninteracting particles [12], assume a simple form: they are pt{l) times 
the difference between the left and right corresponding injection rates, e, 7, respectively, 
namely 

Jp = PtiTL-TR) : Ju = A(eL-eR) • (9) 
Using (8) to eliminate 7 in favor of /i we obtain, 

_ Xpt{l)l (j{D+l)/2 a/T\ 

{Inmyi^'A^ J' ^^^^ 



Taking total differentials of (10) in the variables XjT and /i/T and comparing the re- 
sulting expression with Eq. (4) we obtain exact microscopic expressions for the Onsager 



coefficients and thus, for the TE transport coefficients, namely 



(7 = 



5 = - 



K = 




PT 



1/2 



(11) 



Note that for a chaotic billiard channel with a diffusive dynamics, the transmission 
probability decays as /??(/) °^ which means that all the elements of the Onsager 
matrix L become size independent. 

Finally, plugging (11) into (1) and noting that Cy = D/2 is the dimensionless heat 
capacity at constant volume of the gas, we obtain 



where for simplicity we have called cy = Cy + 1/2. A particular case of (12) was 
previously obtained, for noninteracting monoatomic ideal gases in 3 dimensions [13]. 

The analytical expressions for the TE transport (11) and (12) are exact [1]. They pre- 
dict and increase of the TE efficiency with D. To verify this and further compare in 
the next section with the interacting particles, we have considered a gas of composite 
particles with di^t > 1 internal rotational d.o.f.[l~\ inside a Lorentz gas channel (Fig. 1) 
coupled to two thermochemical baths. The model of composite particles and its dynam- 
ics is explained in the appendix B. We have numerically measured the TE efficiency t] 
as the ratio 



where A/i = Hr — Hl is the chemical potential difference. For a fixed temperature 
gradient, in Fig. 2 we show 7] as a function of V(/i/r) for different values of D. Note the 
excellent agreement between the numerical data (symbols) and the analytical solution 
of 7], obtained from substituting (10) into (13). In each case, the efficiency reaches 
a maximum value rjmax for some optimal value of V{ii/T). Moreover, rj = occurs 
when Jp = 0, which in turns is determined by the Seebeck coefficient. Consistently with 
Eq. (11), the value of V(/i/r) for which 7] = grows linearly with D. Finally, in the 
inset of Fig. 2 we show the numerically obtained rjmax as a function of D (stars) and 
compare them with the analytical expression of (2), with an excellent agreement. 



We now turn our attention to the interacting case. In absence of particles' interaction, 
ZT is independent of the sample size / and depends on the temperature only through 
the chemical potential term. This is due to the fact that with no interactions, pt depends 
on the geometry of the billiard only. From a physical point of view this means that the 




(12) 



3. PARTICLE-PARTICLE INTERACTION 
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FIGURE 2. Thermodynamic efficiency T] as a function of V [jl/T) for a chain of 15 cells, = 95, 
Tr = 105, and D = 3 (circles), D ~ 7 (squares), D — 11 (diamonds), and D = 15 (triangles). The 
dashed lines correspond to the analytical solution obtained from substituting (10) into (13). In the inset, 
maximum efficiency obtained for the different values of D (crosses), in units of the T7caraot. The dashed 
line corresponds to expression (2). 



mean free path of the gas particles is energy independent. When the particles interacting, 
in general pt depends on the local density and temperature of the gas, leading to a 
more realistic situation [14]. Therefore, the analytical expressions for the TE transport 
coefficients (11) are not longer exact. 

To consider a gas of interacting particles we assume that the lattice disks of mass M, 
freely rotate, so that when a particle collides with a disk, their energies mixes, according 
to the collision rules (26). This induces a particle-particle interaction that is mediated by 
the collision with the lattice disks. The strength of the interaction is determined by the 
mass ratio between disks and particles A = M/m [7]. The noninteracting case is recover 
in the limit of A ^ 0, when the disks are much lighter than the particles. Therefore, we 
can study the effect of the interaction by tuning the value of A. First we have studied 
the TE transport for a gas of monoatomic particles {D = 2). In Fig. 3, we show the TE 
figure-of-merit ZT as a function of A. When A ^ 0, ZT ^3/2, which is the expected 
value of the noninteracting case (see Eq. 12). In the inset of Fig. 3, we also show the 
behaviour of the TE transport coefficients. 

Therefore, Eq. (12) is an upper limit of the interacting case when the interaction 
strength vanishes. This is behaviour is expected to be generic since, as explained in [1], 
any interaction tends to correlate the energy carried by the particles with the external 
gradients, thus decreasing the efficiency. 

Next, we have studied the behaviour of the TE efficiency on D, for and gas of 
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FIGURE 3. Figure of merit ZT (squares) as a function of the mass ratio A for a monoatomic ideal gas 
(D = 2). The Onsager coefficients were obtained from two complementary experiments on a channel of 
L = 10 cells, scatterers of radius R=\, and with parameters: Exp 1) Tl — 900, Tr — 1110, /Xl/7l = 
Mr/7r — —6.3, and Exp 2) 7l = Tr = 1000, ^J.l/Tl ~ —6.1 and ^r/7r = —6.5. In the inset, we show the 
transport coefficients a (circles), K (squares), and (triangles). For the sake of presentation, a has been 
scaled by a factor of 2000. The lines correspond to the scalings A^'/^ (solid), A^'/^ (dashed), and A^^/^ 
(dotted). 



interacting composite particles. In Fig. 4, we show the TE efficiency 7] as a function 
of V(/i/r) for A = 1, for different values of D. As for the noninteracting case, the 
efficiency reaches a maximum value Tjmax for some optimal value of the chemical 
potential gradient. However, consistently with the results of Fig. 3, rjmax of the gas of 
interacting composite particles is smaller than the noninteracting gas. This can be seen 
for small D, in the inset of Fig. 4 were rj^ax is plotted for the considered values of 
D. For more complex particles (larger D), rjmax approaches the noninteracting solution. 
Although the behaviour at large D is particular of the rotating-disk interaction, we cannot 
discard that the same could be observed for a more generic type of interactions. Indeed, 
when D ^ I, the amount of energy of the lattice disk becomes negligible compare 
to the total energy contained in the composite particle. Furthermore, for the gas of 
interacting particles, rj also grows with D, indicating that the microscopic mechanism 
for the increase of the TE efficiency discovered in [1] also applies to the interacting case. 

4. CONCLUSIONS 

We have studied the effect of interaction on the properties of thermoelectric transport 
of gases of polyatomic molecules. We have shown that while the TE efficiency of gases 
of noninteracting molecules is an upper bound for the efficiency of the interacting case 
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FIGURE 4. Thermodynamic efficiency T] as a function of V (/x/T) for a chain of 15 cells with A = 1, 
Ti — 95, Tr = 105, and D = 3 (circles), D = 1 (squares), D = 11 (diamonds), and D = 15 (triangles). 
The dashed lines correspond to the analytical solution obtained from substituting (10) into (13). In the 
inset, maximum efficiency obtained for the different values of D (crosses), in units of the TJcarnot- Here, the 
results for D = 2,3 and 4 are also shown. The dashed line corresponds to the analytical expression (2) for 
a gas of noninteracting particles. 



when the interaction strength vanishes, the later still increases with the complexity of 
the molecules. 



A. COUPLING WITH THE THERMO-CHEMICAL BATHS 

Consider a infinite 2-dimensional particle reservoir at equilibrium temperature T and 
density p, filled with point-like composite particles of mass m and dint d.o.f.. The 
reservoir, coupled to the system through an opening of section A, exchanges particles 
with the system, so that in the neighbourhood of the coupling the system is at equilibrium 
with the reservoir. Our aim is to obtain the rates at which particles and energy are 
injected into the system. Inside the reservoir each component of the velocity of the 
particles is distributed according to the Maxwell-Boltzmann distribution 

fr{v)dv = 2n(^:^) v^^-'e^^dv. (14) 

Note that (14) is independent of dint- To obtain the particle injection rate 7 we need to 
compute how many particles of the reservoir hit the opening per unit time. From all the 
particles moving in a given direction v, the number of particles that cross the opening in 
an infinitesimal time interval dt, are those contained in the parallelepid of cross section 
A and length vco^{d)dt, namely pXvco^{9)dt. 



Since the reservoir is at equilibrium, the probability that a particle with speed between 
V and v + dvis injected into the system in a time interval dt is obtained as 



W{v)dvdt = / # / dOv^sinO {pXv cose) {2k) ^ fT{v)dvdt , 

Jo Jo (15) 

Taking v = ^/lE/m in (15), where £ is the translational energy, one obtains the proba- 
bility that a particle with energy E between E and E + dE is injected into the system in 
a time interval dt as 

In terms of W{E), the particle injection rate is defined as 7 = jQW(E)dE, yielding 

The energy injection rate is obtained as £ = /(£'), where {E) is the mean energy of the 
injected particles, namely {E) = (Etrans) + {Eim), where (Etrms) is the mean energy of 
the translational d.o.f. given by 

JlEW{E)dEdl_3 

{t^trans)- g^^E)dEdt ~ 2^ ' 

and {Eint) = ^kT is the mean energy of the internal d.o.f. Note that the mean transla- 
tional energy of the injected particles is not ^kT but ^^kT. Denoting the total number 
of d.o.f. as D = d + dint we finally obtain 

^ p{kTfl\ (18) 
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B. COLLISION RULES 



In this appendix we derive the collision rules that mix the energy among the translational 
components of the particle's velocity, its internal degrees of freedom. We consider 
composite particles with d^t internal rotational J.o./Each "particle" of mass m can be 
imagined as a stack of Jint small identical disks of mass m / Jint and radius r<^R, rotating 
freely and independently at a constant angular velocity (Oi, i— I,.. .,dint. The center of 
mass of the particle moves with velocity v = (vj^, v^). 

We assume that at a collision of a particle with a lattice disk, the normal component of 
the particle's velocity is reflected, and the tangent component, and internal angular ve- 
locities change, so that the total energy and local momentum are conserved. Introducing 



the following notation 

= a^'^coi , for 1=1,..., dint , 

where oc = 0/m and = mr^ /Idint is the moment of inertia of the internal disks, the 
particle's energy becomes 

E='^{^l + t^f\ ■ (20) 



Denoting with primed (unprimed) variables the velocities after (before) the collision, the 
equation for the energy conservation is 

E^'? = E^f, (21) 

i=0 (=0 
and for the conservation of local momentum 

dint difit 

^(, + a'/^l^^'i = ^o + a'/'l^^i. (22) 

i=i i=i 

To solve for the primed momenta, we cast (22) as 

^0 = ^0 - dintJ^ and ^/ = ^i + -Lj(r. (23) 
Substituting (23) into (21) we obtain for the coUisional torque ^ 

dint {I + dint a) \ / 
Substituting (24) back into (23) and using (19) the we obtain the collision rules 



\-dinta la 



l^dinto'^ ' l+dintaf^^^' (25) 
, la 2 ^ 

CO: = COi + ; Vt ; — Tz ; r ) 6), . 

' ^ l+dinta dint{l+dinta)^^ 

For the interacting case, namely when the lattice disks rotate, these can be considered 
as one additional rotor of mass M and radius R and angular velocity c3. Following the 



I + dint ah 



same procedure, the reader can easily verify that the collision rules are 




) 



(26) 



where A = M/m and 



^ = 2(4. + ^ + 



1 



) 



(27) 



These collision rules are a generalization of the ones introduced in [7]. Thus, they are 
deterministic, time reversible and preserve the energy and local angular momentum. 



1. G. Casati, C. Mejia-Monasterio, and T. Prosen, Phys. Rev. Lett. 101, 016601 (2008). 

2. G. Mahan, B. Sales, J. Sharp, Phys. Today 50, 42 (March 1997). 

3. A. Majumdar, Science 303, 778 ( 2004). 

4. M. S. Dresselhaus et al. Adv. Mater 19, 1043-1053 (2007). 

5. T. E. Humphrey and H. Linke, Phys. Rev. Lett. 94, 096601 (2005). 

6. C. Maes and M. H. van Wieren, J. Phys. A: Math. Gen. 38(2005) 1005; J. Volhner, T. Tell, and L. 
Matyas, J. Stat. Phys. 101 (2000) 79; C. Van den Broeck, in Advances in Chemical Physics Vol.135 
Eds. S. A. Rice, (2007 John Wiley & Sons, Inc.). 

7. C. Mejia-Monasterio, H. Larralde, and F. Leyvraz, Phys. Rev. Lett. 86 (2001) 5417. 

8. D. J. Bergman and O. Levy, Appl. Phys. 70 (1991) 6821. 

9. R L. Walstrom, Am. J. Phys. 56 (1988), 890. See also H. Larralde, F. Leyvraz, and C. Mejia- 
Monasterio, J. Stat. Phys. 113 (2003), 197, where the equivalence between n and He has been studied 

for a similar model to the one we study here. 

10. C. A. Domenicah, Rev. Mod. Phys. 26 (1954), 237. 

11. H. B. Callen, Phys. Rev. 73 (1948) 1349; S. R. de Groot and R Mazur, Non-equilibrium Thermody- 
namics (Dover, New York, 1984). 

12. Assuming that particles do not interact means that we fully neglect the phonon transport. 

13. C. B. Vining, Mat. Res. Soc. Symp. Proc. 478 (1997) 3. 

14. J.-R Eckmann, C. Mejia-Monasterio, and E. Zabey, J. Stat. Phys. 123 (2006) 1339. 



REFERENCES 



